Universality class of holographic superconductors 



O 



Kongo 

Department of Engineering, Shibaura Institute of Technology, Saitama, 330-8570, Japan 

Makoto Natsuuni(|l| 

Theory Division, Institute of Particle and Nuclear Studies, 
KEK, High Energy Accelerator Research Organization, Tsukuba, Ibaraki, 305-0801, Japan 

Takashi OkamurgH 

Department of Physics, Kwansei Gakuin University, Sanda, Hyogo, 669-1337, Japan 

(Dated: May 17, 2010) 

We study "holographic superconductors" in various spacetime dimensions. We compute most 
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I. INTRODUCTION 

The AdS/CFT duality [H-ll] has been a useful tool 
to get insights for QCD. (See Ref. @ for a review.) 
Recently, there have been many attempts to apply the 
AdS/CFT duality to condensed- matter physics. (See 
Ref. for a review.) The critical phenomena may 
be useful to make further progress in such applications. 
Namely, given a field theory system, the universality class 
of the system may be helpful to identify the dual gravity 
system. In this paper, wc study both the static universal- 
ity class and the dynamic universality class of holographic 
superconductors in various dimensions.^ 

A holographic superconductor is a solution of an 
Einstein-Maxwell-scalar system [9l-[l]| or an Einstein- 
Yang-Mills system [l^]. There are two branches of the 
solution. For T > Tc, the solution is the standard 
{p -\- 2)-dimcnsional Reissner-Nordstrom-AdS black hole 
(RN-AdSp+2) with scalar * = 0. For T < T^, the so- 
lution becomes unstable and is replaced by a charged 
black hole with a scalar "hair." The analytic form of the 
solution is unknown though, so one often uses numerical 
computations or uses the probe approximation, where 
the backreaction of the scalar field on the metric can be 
neglected. 

The analysis of these papers indicates that a holo- 
graphic superconductor has the conventional mean-field 
behavior, but one had better compute all critical expo- 
nents of the system to confirm this. This is our main 
purpose. One might argue that it is not necessary to 
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compute all critical exponents since they arc related by 
scaling relations. However, it is not a priori obvious that 
a holographic superconductor obeys the scaling relations 
since it is not a standard statistical system but is a grav- 
ity system. Thus, computing these critical exponents is 
also important in order to check the AdS / CFT duality. 

In this paper, we focus on the high-temperature phase, 
which simplifies the analysis. We study the linear per- 
turbations of the bulk equations of motion. Our main 
results are summarized as follows: 

1. The static universality class of holographic super- 
conductors is the one for the Gaussian fixed point. 
In other words, the static critical exponents take 
the standard mean-field values. 

2. The dynamic universality class is model A with dy- 
namic critical exponent z = 2 satisfying the dy- 
namic scaling relation z = 2 — 77. 

The results are obtained by analytically studying the 
equations of motion without the probe approximation 
(Sec. IIV[) . These results are also confirmed in numerical 
computations with the probe approximation (Sec. \V\ . 

In the static critical phenomena, there are six critical 
exponents (a, /?, 7, (5, i^, 77). (See App. HIl) Wc compute 
these exponents except /3 and d. (The exponent /? ap- 
pears only in the low-temperature phase, and 5 is beyond 
the scope of the linear perturbations. See Sec. [IV CI ) The 
critical exponents we found are natural since fluctuations 
are suppressed at large- so that the mean-field values 
become exact. Moreover, our results are independent of 
spatial dimensionality which is typical for mean-field re- 
sults. 

In the dynamic case, the effect of diverging correla- 
tion length ^ also appears in the relaxation timc^ and 
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in the transport coefficients, wfiicfi is known as tlie criti- 
cal slowing down. The divergence is parametrized by the 
dynamic critical exponent z. 

The dynamic critical phenomena in the AdS / CFT du- 
ality has been studied in Ref. [l^l using single R-charge 
black holes in various dimensions |2ll - l23l | . The order pa- 
rameter of the transition is the charge density p. Since 
it is a conserved quantity, the simplest possibility for the 
dynamic universality class is model B according to the 
classification by Hohenberg and Halperin [23]i and in- 
deed it has been shown to be model B. 

For the holographic superconductors, the order param- 
eter of the phase transition is the scalar condensate (O) 
which is dual to the scalar field . This is not a conserved 
quantity, so the model should belong to model A. 

We set up our notations and conventions in Sec. HIl 
and set up the perturbation equation in Sec. IIIII In 
Sec. IIVI we formally solve the perturbation equation as 
a double-series expansion in (a;,^), where w,g are the 
frequency and the wave number, respectively. We argue 
that, rather generically, the values of these exponents can 
be understood simply in terms of a double-series expan- 
sion, and we sec no sign which indicates the expansion 
breaks down. (If the exponents took the different values, 
this should appear as the break-down of the expansion 
somewhere.) Our argument is rather generic so that it 
may cover other, yet to be discovered, systems with the 
mean-field behavior. However, the argument in Sec. IIVI 
relies on certain assumptions, so in Sec.|V]we numerically 
solve the perturbation equation in the probe approxi- 
mation and obtain the exponents for the special case of 
the RN-AdSs solution with scalar mass I'^rv? = —3. In 
App. [XI we review the basics of critical phenomena. In 
App. [El we show the anomalous dimension 77 = using 
a similar technique as in Sec. IIVI 



The equations of motion are given by 



^ 2A - — + y 



= V^i^^'^ - f = 



= Di'B,,^ f-r^* 



where the current j is defined by 



(2.5a) 

-gF^-) - 3^ , (2.5b) 
(2.5c) 



--^^D^^) . (2.6) 



A. Thermal equilibrium solution 

We study a static asymptotically AdS black hole with 
planar horizon, and we take the following ansatz: 

dsi^^ = — ( : de C H— (u) dxp^ 



/(") 

A^, = -*(u) (dt)^, , 



(2.7a) 

(2.7b) 
(2.7c) 



where C is a constant which is related to the black hole 
temperature [See Eq. (|2.15|) ] . Without loss of generality, 
one can choose 



f{u = 0) = H{u = 0) = 1 



(2.8) 



II. HOLOGRAPHIC SUPERCONDUCTORS 

In order to study the critical phenomena of holographic 
superconductors, we consider a (p + 2)-dimensional 
Einstein-Maxwcll-scalar system: 



C 



R~2A 



vm), (2.1) 



where the covariant derivative D^, a cosmological con- 
stant A, and the scalar potential Fd^*!) are given by 



2A = - 



ieA,, 



p{p+l) 



(2.2) 

(2.3) 
(2.4) 



peared in the second-order hydrodynamics Il3ll . The AdS/CFT 
duality for the second-order hydrodynamics was developed by 
Refs. [l3-[Ii. 



This w = is the AdS boundary. The horizon is a so- 
lution of f{u) = 0, and we can set that the horizon is 
located at u = 1 by an appropriate scaling. 

The ansatz is partly motivated by R-chargcd black 
holes [2T| - [23j . When p = 2,3, the ansatz corresponds 
to the RN-AdSp+2 limit of R-charged black holes. (The 
p = 2 and p = 3 R-charged black holes have four U{1) 
and three U{1) charges respectively, and the RN-AdS 
limit corresponds to the equal-charge cases.) 

The equations of motion reduce to 



= (Einstein equation) , 

^ d f d Pm'^H^ _ e^H'^ ^2 



0= 



du uP du uP+'^ 



CuPf 



(2.9a) 



(2.9b) 



0=|^:AiA_?!!£!^|^|^U, (2.9c) 



H"^ du uP du u^ f 

^ d^\ ^.d* 

du du 



(2.9d) 
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Equation (j2.9dp implies that the phase of ^ must be 
constant so that one can set ^ to be real without loss of 
generality. 

The background solution is obtained from Eqs. (|2.9p 
by imposing (i) the regularity condition at the horizon 
and (ii) the asymptotically AdS condition. The former 
is given by 



f{u = 1) = ^{u = 1) = , 
H{u = 1), \I>(?i = 1) = const. , 

and the latter is given by Eq. (|2.8|) and 
where A± are defined by 




(2.10a) 
(2.10b) 



(2.11a) 
(2.11b) 



(2.12) 



Here, fj, is the R-charge chemical potential and C corre- 
sponds to the R-charge density. 

References [l3, EH have solved Eqs. (|2.9I) by imposing 
the boundary condition -i/j^"-* = at the AdS boundary, 
and they found that the solution has two "phases" : (i) a 
black hole solution with no scalar hair at high tempera- 
tures and (ii) a black hole solution with a scalar hair at 
low temperatures. According to the standard AdS/CFT 
dictionary, -i/j'-'''-' corresponds to the expectation value of 
the scalar operator dual to the scalar field, while ■j/)'-"-' 
corresponds to the dual source.'^ This suggests that the 
hairy solution with i/)'^' corresponds to a sponta- 
neous condensation of the dual operator. Furthermore, 
the hairy solution has the expected behavior for a super- 
conducting phase, i.e., (i) the divergence of the R-charge 
(DC) conductivity and (ii) an energy gap proportional to 
the size of the condensate. 



B. High-temperature phase 

At high temperatures, the black hole has no scalar 
hair ('S' = 0), and the solution is given by the RN-AdS 
solution: 



H = l + KuP-^ , 

f - {1 + k)^uP+^ , 

1 + K 



At = /i 1 



H 



(2.13a) 
(2.13b) 

(2.13c) 



This background (|2.7|) with Eqs. (|2.13p is parametrized 
by three parameters k, ^, and (. The RN-AdS black 
hole may be parametrized by the temperature T and the 
chemical potential fi (in the grand canonical ensemble). 
Thus, only two among three parameters are independent. 
In fact, they are related by 



^ = lC\h^^ k'/'(1 + k)^ . (2.14) 

V 1 

Clearly, the Schwarzschild-AdS solution (SAdS) corre- 
sponds to K = 0. The black hole temperature is related 
to C: 

2.T = C(l + .)^ P+1~(P~1)" . (2.15) 



In particular, when k = 0, C has a simple relation with 
the SAdS temperature Tq: 



c 



p + l 



-TrTn 



(2.16) 



The charge density p of the RN-AdS black hole is given 

by 



lp={p-l){l + K){K)P-'tJi 



(2.17) 



III. PERTURBATION EQUATION AT HIGH T 

Our aim is to investigate both the static and the dy- 
namic critical phenomena for holographic superconduc- 
tors. The basic information about these phenomena such 
as critical exponents can be read off from the two-point 
correlation function. In the AdS/CFT duality, the cor- 
relation function of operator O in the boundary theory 
can be obtained from the bulk perturbation of the dual 
scalar field. 

We study the critical phenomena of holographic super- 
conductors approaching from high temperature. This has 
two advantages: (i) The background solution is analyti- 
cally known; (ii) The scalar perturbation -0 := 4' — '4' de- 
couples from the electromagnetic and gravitational per- 
turbations because the background solution has no scalar 
hair. 

Under the ansatz ip = V'ro,q(")e~''^*^*'^'^, 0ro,q(u) obeys 



0=(uP±Ld 
\ du uP du 



2p 

H~ 9 

(tt3 + a)' 



/ 



V'ro,q(u) , (3.1) 



When -(p -I- 1)^/4 < Z^m^ < -{p + 1)^/4 -|- 1, both operators 
are normahzable [Ssll so that we have two choices of the order 
parameter, (Oi) and {O2) in the notations of Refs. [Tol. In 
those cases, we consider only the operator {02)- 



where we introduced the following dimensionless quanti- 
ties: 



oj (p + 1)lu 
ir := — = — — 

C 4^ro 



\q\ {p+l)\q\ 

1 •= ^ ^ A rr ' 3.2) 

C 47rTo 
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and 



eAt ( 1 
C Vl + '^ 

<t:=(1 + ^)^. 



(3.3) 
(3.4) 



Roughly speaking, a cx /i/T. Our background may be 
written by T and /i, so cr is a unique dimensionlcss pa- 
rameter to characterize the background. Also, note that 
a and k are related by Eq. p.l4p : 



yi^(/e).i/2(i + ,)-^ . (3.5) 



Equation p.ip is characterized by six dimensionlcss 
parameters: 

• three parameters p, (^to), and (Ze) which 
parametrize the action, 

• a parameter cr (or k) which characterizes the back- 
ground, 

• two parameters to and q which parametrize the per- 
turbation. 

The probe approximation often simplifies the analysis. 
In our formulation, the probe approximation is obtained 
by taking e — > oo keeping a fixed. Equation p.Sp tells 
that K cx (le)~'^ ^> in this limit, namely the background 
(|2.7p becomes the SAdS solution. This limit may superfi- 
cially look different from the conventional one employed 
in Ref. [13, where one takes e — oo keeping (e^) and 
(e$) fixed. But both of these limits yield the SAdS so- 
lution. 

In this approximation, Eq. (j3.ip is characterized by five 
parameters: p,{lm),a (with k = 0), Id, and q. We will 
employ the probe approximation for numerical analysis 
in Sec. |Vl but we will not employ the approximation for 
the analytic argument in Sec. IIVI 

Near the horizon w ~ 1, Eq. p.ip becomes 

" ^ (^ll^ll + + ^)^lr2) r^r^A^) , (3.6) 

so the solution is given by ^pm.q ^ (1 — u)^'tIt. We 
impose the "incoming wave" boundary condition at the 
horizon, which corresponds to the retarded condition: 

^ro,q ^ (1 - u)-'^ . (3.7) 

Near the AdS boundary u ^ 0, Eq. p.ip becomes 



= ( uP—u-P— - 
ou ou 



(3.8) 



and ^AtD.q has the same fall-off behavior as the background 



*: 



(3.9) 



The critical phenomena near the second-order phase 
transition can be extracted from the two-point correla- 
tion function or the response function. In the AdS/CFT 

duality, the source term corresponds to V'ro.q: ^-i^d the 
dual order parameter expectation value {Oxn,q) corre- 
sponds to V'rotq J so the response function is given by 



tro,q 



yro,q 
rro,q 



(3.10) 



Thus, our task is to solve Eq. p.ip under the boundary 
condition p.Tp . obtain the coefficients i/'ro^q in Eq- p.9p . 



and study the behavior of the response function Xm,q 
4^1 1 



IV. GENERAL STRUCTURE OF RESPONSE 
FUNCTION BY (ro, q)-EXPANSION 



First, we formally solve the perturbation equation p.ip 
as a double-series expansion in (tri,q). We demonstrate 
that the critical exponents for holographic superconduc- 
tors can be understood simply in terms of a double-series 
expansion. 

In order to implement the "incoming wave" boundary 
condition p.7p . it is convenient to introduce a new vari- 
able Lp as 



'0ro,q(w) =: T{u)ipm,(\{u) 



(4.1) 



The choice of the function T{u) is arbitrary as long as it 
satisfies 



but we choose 



for u ~ 1 



T = exp 



I I du — - — (ir + 21) 

/ 



(4.2) 



(4.3) 



Then, the "incoming wave" boundary condition at the 
horizon corresponds to (^,^_q(u = 1) = const., and 
Eq. (|3.ip becomes 







+ Bi{u)— + So(w) Vro,q(w) 1 



d f H~ 
Sl(w) = £ln^ + 2^^(m+2l) 



Near the AdS boundary, ipm.q behaves as 

/ \ (-) A_ , (+) A+ 

¥'ro,q(w) ^ <(5ro,q M + fm/q U + 



(4.4a) 
(4.4b) 



(4.4c) 



(4.5) 
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where A± (A_ < A+) is defined in Eq. ((2?T2l) . 

Let us solve Eq. (|4.4ap as a double-series expansion in 
(tt',q): 

+ 0{w\K>q\q^) . (4.6) 
Each variable obeys the following differential equations: 



Cipo{u) :-- 



^<P(i,o)(w) = -^(l,0)(") 



ipo{u) = , 

(4.7a) 
(4.7b) 
(4.7c) 



where Ci(u), Co(u), J(io)(m), and J(-oi)(u) are defined 

by 



d f H~ 
Ci(zi) = -f ln4 + 2z— ^21 

au uP J 



(4.8a) 



Co(u) = -?W:^+z^A (^2t) , (4.8b) 



u^/ f du \ uP 



and 



-^(1,0) (w) = 



d d I H— 
2— + — ln- 



du du \ uP 



(4.9a) 
(4.9b) 



because the thermodynamic susceptibility x is the re- 
sponse function p.lOp for the stationary homogeneous 
source x = Xro=ci=o- 

In the high-temperature phase T > Tc, it is possible to 

show that the order parameter '^q'*'' vanishes if there is 
no source (^q Thus, </jg^'' is proportional to (/?q and 
X is finite, which is the behavior similar to paramagnets. 

On the other hand, at the critical point, the sponta- 
neous condensation occurs, so one must have (/Jq"-* It,, = 

and <y3o^''|Tc 7^ 0; we assume that this is the case. This 
is nothing but the divergence of the thermodynamic sus- 
ceptibility X at the critical point. Then, the critical ex- 
ponent of the thermodynamic susceptibility can be ex- 
tracted from the behavior of (^g As T approaches 
from the above, we assume 



lim ibr, 

T\Ta 



(-) 







lim i/ir. 



(+) 



finite . 



(4.13) 



So, denote the deviation from the critical temperature Tc 
by ct := T/Tc - 1 > 0, and write i/ig"-* as 

4'^ = cx (7 > 0) . (4.14) 

Then, the thermodynamic susceptibility behaves as 

X (X ^^+V^[)"^ « (4-15) 



namely 7 in Eq. (|4.14p is indeed the critical exponent of 
the thermodynamic susceptibility. 



B. First order in (ro, q)-expansion 



A. Zeroth order in (ro, q)-expansion 

From the boundary condition, the zeroth-order solu- 
tion ipoiu) must be regular at the horizon, and it must 
satisfy 



u=l 



Co 
Ci 



u=l 



/' 



u=l 

(4.10) 



This determines the boundary condition for (/Jq up to an 
overall constant, so we shall set lpq{u = 1) = 1. 

Suppose that one solves the differential equation (|4.7ap 
under this boundary condition by integrating out from 
the horizon to the AdS boundary. Then, one can deter- 
mine the behavior of ^po{u) as 



(4.11) 



and determine the ratio (ySg^^/t^g \ This ratio is propor- 
tional to the thermodynamic susceptibility of the order 
parameter, x'- 



(4.12) 



It is easy to solve Eqs. (j4.7bp and (|4.7cp formally. Let 
(po be a zeroth-order solution which is linearly indepen- 
dent from (po- Then, one can construct the Green func- 
tion of the operator C which satisfies the regularity con- 
dition at the horizon: 



G{u, u ) := 6{u — u )lpq(u)(Po{u' ) 
+ d{u' - u)(po{u)ipo{u' ) , 



(4.16) 



W{u) := — (fioiu) ~ ipo{u)— , (4.17) 

du du 

where the Green function G and the Wronskian W satisfy 
CuGiu,u) = W{u)S{u-u'), (4.18) 



— W{u) + Ci{u)W{u) = 
du 



(4.19) 



Using the Green function (|4.18p . the solution of 
Eq. (|4.7bp is expressed formally as 

= du G{u,u ) , 

, . r , , ,/(i,o)(^') 

= ipo[u) / du ipo{u )- 



JO Wiu') 
+ Mu) £ du'Mu' )^^^^|^ . (4-20) 
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and similarly for Eq. (|4.7c|) . 

We also assume that the first-order solutions have the 
same fall-off behavior as Eq. ()3.9p . Then, 



(-)„A_ , (+), A + 



Comparing this with Eq. (j4.5p . wc obtain 

'/'tU = v'^^l = "^o"^ + "^'/'[r]) + qV(o,l) ' (4-22a) 
^Ltq - ^i^q = V'^^' + "'^(to) + ■ (4.22b) 



C. The order parameter response function 

For readers' convenience, let us list our assumptions 
used so far: 

1. The first-order solutions V'(i,o) ^^'^ V(o.i) indeed ex- 
ist, 

2. They have the same fall-off behavior as Eq. p.9p . 

3. The singular behaviors at the critical point come 
only from (^g i.e., (^q 'It^ ~ 0. In particular, we 
shaU use ^^o'\T^,v\il)\T^,f\^l)\T^ 0. 

Using the above assumptions and the results in previous 
subsections, the response function p.lOp is given by 

ifiQ + '"V5(l,0) ^ 1 ^(0,1) 



n + tt"/'(i,o) + qV(o;i) 

^(0,1) ' 



(4.23) 



where c and ^ are defined by 



c *^(i,o)/^(o7d ' ^= ^(o,V^o"^ ■ (4-24) 

In particular, the static and dynamic response function 
are given by 



'''(0,1) ^ 



(4.25) 
(4.26) 



From the response functions (|4.23p and (|4.25p . we can 
extract the following information: 



• The form of the static response function (|4.25p tells 
that ^ is indeed the correlation length, and ^ di- 
verges as ^ oc e^^^^ from Eqs. (14. 14^ and (|4.24l) . 



This determines the critical exponent v (defined by 
^ cx e^"^) as = 7/2. 

At the critical point (ey = 0), the static response 
function behaves as XrD=o,q|Tc q"^. This means 
that the anomalous dimension 77, which is defined 
by XrD=o,q|Tc q''"^, vanishes. These critical expo- 
nents 7, V, and 77 satisfy the static scaling relation 
7 = v{2 — 77), as expected. 



• The dynamic response function (j4.26p indicates 
that the homogeneous perturbation of the dual or- 
der parameter decays on the time scale rq=o ^ c^^^ 
This is the critical slowing down, in which the re- 
laxation time diverges as a power of the correlation 
length (if we assume c 7^ near the critical point). 
The dynamic critical exponent z, which is defined 
by rq=o c)C is equal to 2. This exponent z also 
appears in the dynamic response function at the 
critical point: it appears in the q-dependence as 
Tqlr^ cx q~^. In fact, Eq. (|4.23p gives ^ cq"^ at 
the critical point (when ^ = 00). 

This is our main result. 

The explicit value of the exponent 7 is determined ei- 
ther from the numerical computation in Sec. |V] or from 
an analytic argument in App. [BJ wc find 7 = 1. As a 
consequence, the static critical exponents are 



7=1 



= 7/2 = 1/2 



(4.27) 



which satisfy the static scaling relation 7 — i>{2 — rj). 
Note that our results are independent of spatial di- 
mensionality which is typical for mean-field results. In 
Refs. [2^,1231, the exponent v has been computed in the 
low-temperature phase. We find that the exponent takes 
the same value in the high-temperature phase, which is 
typical in critical phenomena. 

It is important to check the other static scaling rela- 
tions 



2/3 + 7 = 2 



7 = /3(<5-l) 



(4.28) 



Here, a is the exponent for the heat capacity, S is the ex- 
ponent which gives the behavior of the order parameter 
at the critical point: ( Orn=q=o ) oc 

IV'^V^''- The expo- 
nent /3 determines the behavior of the order parameter 
in the ordered phase: (01^=^=0) IctI^- We focus on 
the high-temperature phase, so we cannot obtain f3, but 
^ = 1/2 from Refs. [iMl- 

As we approach the critical point from the high- 
temperature phase, the bulk spacetime is the RN-AdS, 
so the heat capacity is given by the one for the RN-AdS. 
Thus, the heat capacity converges to a constant value 
as we approach the critical point, which implies a = 0. 
Hence, the first scaling relation in Eq. (|4.28|) is confirmed. 

In order to check the second scaling relation in 
Eq. ()4.28|) . one needs to compute S, but this analysis 
itself does not give S. The would-be value for 5 is (5 = 3, 
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and this means that V^g"'' c>c (ipo^^)^, so the hnear pertur- 
bation is not enough to derive the result. (See however 
Ref. , where the authors show that the standard Lan- 
dau potential fits their numerical results. This indicates 
(5 = 3.) 

We find that the static critical exponents take the stan- 
dard mean-field values (except 6). One might suspect 
that our argument is simply a consequence of double- 
series expansion. But the whole point is that the uni- 
versality class of holographic superconductors can be un- 
derstood in this way. If the exponents took the different 
values, this should appear as the breakdown of the ex- 
pansion somewhere, but we see no sign of the breakdown. 
However, our argument relies on certain assumptions, so 
we numerically obtain the exponents in the next section. 

Now, let us turn to the dynamic critical exponent z. 
We consider the linear perturbations of the bulk field. In 
the dual field theory, this corresponds to the following 
approximations: 

(a) One ignores the couplings between the fluctuations 
of the order parameter with different wave number 
(the Gaussian approximation for the fluctuations). 

(b) One ignores the nonlinear couplings between hydro- 
dynamic modes and the fluctuations of the order 
parameter (no mode coupling). 

This approximation is known as van Hove theory ("con- 
ventional theory") of the dynamic critical phenomena, 
which predicts that the dynamic critical exponent z is 
related to the static critical exponent 77 (anomalous di- 
mension) hy z = 2 + a — rj. Here, a = when the order 
parameter is not a conserved quantity (model A) and 
a = 2 when it is a conserved one (model B). 

In holographic superconductors, the order parameter 
is the scalar condensate (O) which is not a conserved 
quantity, so we expect z — 2 ~ j] from van Hove theory. 
This is consistent with our results 77 = 0, z = 2 obtained 
from the (tr, q)-expansion, implying that the AdS/CFT 
duality also holds in the dynamic case. In other words, 
these black holes do obey the theory of dynamic critical 
phenomena. 

It is interesting to see whether z obtained from nonlin- 
ear perturbations is consistent with the properties pre- 
dicted by the dynamic renormalization group analysis. 
Under nonlinear perturbations, one generally expects (a') 
the violation of the Gaussian approximation and (b') the 
appearance of mode-mode coupling. Since the dual field 
theory is at large- A^, the effect of renormalization from 
short-wavelength fluctuations would be suppressed and 
we could ignore (a'). On the other hand, (b') has a sub- 
stantial effect, which would change the dynamic univer- 
sality class from model A (maybe to model F which is 
the universality class of '''He) . It would be interesting if 
one could obtain z in the gravity side which is consistent 
with the dynamic renormalization group analysis. This 
could become a strong support of the AdS/CFT duality 



in the dynamic regime.' 



V. NUMERICAL RESULTS 

The results in Sec. llVl and App. [Blare obtained by as- 
suming that (i) Eq. (j4.4p can be solved as a double-series 
expansion in (tn, q), (ii) each other solution preserves the 
fall-off behavior p.9p . and (iii) singular behaviors come 
only from (p^ \ In this section, we solve Eq. (j4.4[) numer- 
ically for a particular set of parameters and show that 
our results based on the above assumptions are indeed 
correct. 

We numerically solve Eq. ()4.4p under the boundary 
condition 



^ 1 



Bo 
Bi 



(5.1) 



.(+) /,„(-) 



from 



and obtain the response function Xm,c\ oc 'p)x,,ql 'rm,q 
the asymptotic behavior (|4.5p . 

To simplify the analysis, we perform the numerical cal- 
culation in the probe approximation. The perturbed so- 
lution then depends on flve parameters p, (Im), cr, tr, 
and q (Sec. IIII[) . As an explicit example, we consider 
five-dimensional bulk spacetime {p = 3) with the scalar 



mass Pm^ = —3. which was discussed in Ref. [27[. In the 



probe approximation, the background (|2.7p becomes the 
SAdS solution, so/ = l-u'*,i? = l. 



A. Critical point and thermodynamic susceptibility 

In our formulation, the dimensionless parameter a de- 
termines the phase structure, and let us first find the 
critical point CTc. We solve 

• EOM: Eq. (|Oa|) with ir = q = 0, or Eq. dUTaj), 

• Boundary conditions: Eq. (jS.ip at the horizon and 
(^Q ^ = (corresponding to a spontaneous conden- 
sate) at the AdS boundary. 

Then, (Tc is obtained as an eigenvalue problem under two 
boundary conditions. We use a numerical method devel- 
oped by Horowitz and Hubeny to solve such an eigenvalue 
problem psj . The parameter is numerically found as^ 



(7^^- 4.15686 or T^- 0.19797 x 



ep\ 1/3 
21- 



(5.2) 



* This nonlinear effect (b') may be subleading in the 1/A'^- 
cxpansion. We thank Misha Stephanov for pointing this out 
to us. 

5 Reference ^ found that Tc = 0.198 X pj/^ in units I = e = 1, 
which coincides with our numerical result. Note that the R- 
chargc density Phr in ^ei. [StI is related to our p by pjjR = p/2. 
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FIG. 1: (color online). The value of ip[ ^ ("+", (red)) and 
(¥){,+' - 4.585)/100 ("x", (blue)) as a function of a. As cr 



approaches the critical value <Jc ~ 4.15686, (ySg ' goes to zero 



and ly'o^' goes to a finite constant. 
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1 
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FIG. 2: (color online). The thermodynamic susceptibility x 
as a function of \~o ja^- Plotted are 1/x for = 10~^n 

(n = 0,l,-- - ,20). 



To examine the critical behavior of the thermodynamic 
susceptibility deviate cr from the critical point. We 
denote the deviation from the critical point by eo-, where 



:= 1 - ojac oc er 
= 10"^ n 



(n = 0,l. 



,20) 



(5.3a) 
(5.3b) 



and see how physical quantities behave under the defor- 
mation. Figure [I] shows <Po^^ as a function of a around 
the critical point (in the high-temperature phase). As 
expected, the critical behavior coincides with Eq. (|4.13p . 
This guarantees that x diverges near the critical point, 
and Fig. [2] shows that x diverges as x oc l/e^^ suggesting 
7 = 1. In fact, the result may be fitted by a polynomial 



(x 10"' 







1 



2 (x10"^) 



FIG. 3: (color online). The correlation length defined 
by ^2 := -1/q*. Plotted are 1/^^ for €^ = W'^n (n = 
0,1,-- - ,20). 



as 

1/x -1.72 X 10"** + 0.25 xea- 4.58 x el 

oc ea , (5.4) 

so 7 = 1 within numerical errors. 

B. Correlation length and static susceptibility 

We then check 77 = 0, = 1/2. These are exponents 
for the static susceptibility Xro=o,q, so it is enough to 
consider the ro = perturbation. 

In order to obtain ly, we solve 

• EOM: Eq. (|44a|) with ir = 0, 

• Boundary conditions: Eq. (j5.ip at the horizon and 



*■ q = (corresponding to no deformation) at 



the Ads boundary. 

This procedure gives the eigenvalue q.^, which corre- 
sponds to the pole of the static susceptibility. The pole 
gives the correlation length ^ by -.^ — 1/q^. Figure [3] 
shows with the interval Aect = 10"^ toward the crit- 

— 1/2 

ical value ctc, which suggests ^ oc ea ■ The result may 
be fitted by a polynomial as 



1/r 1-57 x 10' 

(X , 



13.55 xe„- 11.16 x 



(5.5) 



so J/ = 1/2 within numerical errors. 

The exponent 77 comes from the static susceptibility at 
the critical point XtD=o,q|Tc oc c\^~^, so it is obtained from 
its q-dcpendcnce. We solve 



• EOM: Eq. (|i:ia|) with ir = 0, 

• Boundary condition: Eq. (j5.ip at the horizon. 
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FIG. 4: (color online). The static susceptibility at the critical 
point Xro=o,q\T^- Plotted are l/xit>=o,q It^ for = lO^^n in — 
0,1,-- - ,20). 
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FIG. 5: (color online). The trajectory of the lowest quasinor- 
mal frequency Idqnm for ~ 10~*n (n = 0, 1, • ■ • , 20). As 
a — ;> (Tc, roQNM approaches the origin with equal spacing, and 
it vanishes at a = ac- 



and obtain XiD=o,q|Tc ^fil^lo q/¥'L=o q from the behav- 
ior at the AdS boundary. Figure |4] shows Xm=o,q\T^ as a 
function of q. The result may be fitted by a polynomial 

as 



1/Xro=0,q 



-1.45 X 10"** + 0.018 X q2 + 0.55 x 



cx q 



(5.6) 



so Xro=o,q|Tc oc q ^ within numerical errors. 



C. Relaxation time 

Finally we check z — 2. This is obtained from the low- 
est quasinormal (QN) frequency, which gives the relax- 
ation time. A QN frequency is obtained as an eigenvalue 
by solving 

• EOM: Eq. (|i:ia|) . 

• Boundary conditions: Eq. (|5.ip at the horizon and 

Vn!,l = (corresponding to no deformation) at the 
AdS boundary. 

We expect that holographic superconductors belong to 
model A, so we set q = 0. For the other classes, one in 
general needs to compute Tq in order to obtain z. [As an 
example, see Eq. (|A14p for model B.] Figure [5] shows the 
lowest QN frequency tDqNM on the complex tr-planc with 
the interval Ae^r = 10"^. As a approaches the critical 
value (Tc, the lowest QN frequency approaches the origin 
with equal spacing, and it vanishes at a = ac- The equal 
spacing suggests ttiQNM oc Ca, and in fact the result may 
be fitted by a polynomial as 



«)QNM (0.76- 1.61i) X 10 
- (0.52 + 0.76i) X el 
oc en ■ 



+ (2.23- 0.50i) X 

(5.7) 



The imaginary part of tDqNM is nonpositive, which indi- 
cates that the system is stable in the high-temperature 
phase and is marginally stable at the critical point. Com- 
bining ^ oc with the numerical result (j5.7p . one 

obtains Itiqnm oc or t oc <^^, namely z = 2. 



VI. DISCUSSION 

We find that the critical exponents for holographic su- 
perconductors take the standard mean-field values, but 
it is not a priori obvious. In the AdS/CFT duality, the 
simplest black holes with a second-order phase transition 
are the R-charged black holes. The p = 3 case is dual to 
the J\f = A super- Yang-Mills theory at a finite chemical 
potential, and it has been widely discussed in the litera- 
ture.^ When the black holes have only one R-c harge, t he 
critical exponents take nonstandard values pol . |3i 



(Q;,/3,7,(5)R_chargod 



1 1 1 

2' 2' 2' 



(6.1) 



While this suggests the existence of field theory systems 
with unconventional critical behaviors, their physical in- 
terpretation may not be easy to understand. 

Then, the question is why holographic superconductors 
follow the conventional mean-field behavior, or how one 
can find a gravity dual with the conventional mean-field 
behavior. We do not have a completely satisfying answer, 
but the following argument may suggest a clue. 

Although there are six static critical exponents, only 
two are independent in normal statistical systems due to 



® For example, it has been used to establish the universality of 
the shear viscosity [29l432ll . 
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scaling relations. One would choose a as an independent 
exponent. The heat capacity in many interesting statis- 
tical systems does not have a power-law divergence or 
has only a weak power-law divergence, so a = 0. On the 
other hand, the R-charged black holes have a power-law 
divergence in the heat capacity, and this may be the rea- 
son why the critical exponents for these black holes do 
not take the conventional values. Then, in order to have 
the conventional mean-field behavior, one should prepare 
a gravity system where the gravity sector does not have 
a singular behavior but only a matter sector has a sin- 
gular behavior. A simple way is to pick up a familiar 
black hole solution without a second-order phase tran- 
sition and to couple a matter field just like holographic 
superconductors. If the matter field undergoes a second- 
order phase transition, the whole system is likely to have 
the conventional mean-field behavior. It is clear that the 
gravity sector does not play an essential role in our high- 
temperature analysis or in the probe approximation: it 
simply provides a background. 

From the gravity point of view, one needs a black hole 
which violates the uniqueness theorem (or the no-hair 
theorem) in order to have a second-order phase transi- 
tion. The uniqueness theorem often fails for asymptoti- 
cally AdS black holes, and a simple way to violate it is 
again to add an appropriate matter field. 

The discussion does not ensure that any system which 
satisfies the above conditions has the conventional mean- 
field behavior because we fix only a among two indepen- 
dent critical exponents. But this discussion suggests that 
there may be many more gravity systems with the con- 
ventional mean-field behavior waiting to be discovered. 
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A: Review of critical phenomena 

In this appendix, we briefly review the static and dy- 
namic critical phenomena. For more details, see standard 
textbooks 1351. 



As an example, consider the ferromagnetic phase tran- 
sition. In this case, the magnetization m and the ex- 
ternal magnetic field h are the order parameter and a 
control parameter, respectively. The static critical expo- 
nents (a, /?, 7, (5, i^, ry) are defined by 



Gh (X lerl-" 

m cx \(-t\^ 

XT oc jerr^ 

m cx \h\^l^ 

G(r) cx e-'-/^ 



i oc IetI 



(T < T,) 
(T = T,) 
(T = T,) 



(Al) 



where := {T — T^) /T^^ Ch is the specific heat, xt is 
the magnetic susceptibility, G{r) is the correlation func- 
tion, and ^ is the correlation length, and ds denotes the 
number of spatial dimensions. 

There are six static critical exponents, but not all are 
independent, and they satisfy static scaling relations: 



a -I- 2/3 + 7 = 2 
7 = /3((5~l) 
7 !^(2 - r/) 
2 — a = vd^ 



(A2a) 
(A2b) 
(A2c) 
(A2d) 



Because of these relations, only two are independent 
among six exponents, which suggests that there is some 
structure behind these relations, which is known as the 
scaling law. 

The dynamic universality class depends on additional 
properties of the system which do not affect the static 
universality class. In particular, conservation laws play 
an important role to determine dynamic critical expo- 
nents. As a consequence, even if two systems belong to 
the same static universality class, they may not belong 
to the same dynamic universality class. 

Phcnomenologically, the relaxation to the equilibrium 
is governed by the time-dependent Ginzburg-Landau 
(TDGL) equation. A prototypical example of nonequi- 
librium phenomena is the Brownian motion, which is de- 
scribed by the Langevin equation: 



dv{t) 
dt 



= -Tv{t) + at) 



(A3) 
(A4) 



where F is the friction coefficient, H = /2\s the Hamil- 
tonian, and C,{t) is a random force with {C,{t)) = 0. The 
TDGL equation is the many-body generalization of the 
Langevin equation: 

dyT[\x-y\) +Cit,x) , (A5) 



dt 



Sm{t,y) 



where r(|x — y|) is a transport coefficient which plays the 
role of F, and /[to] is the pscudofrec energy. 
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As an example, consider the Gaussian model with a 
source term: 



I[m; T,h] ~ / dx 



I (Vm)2 



— m 
2 



mh 



(A6) 



where a = ao(T — T^) + • • • and c = cq + ■ ■ ■ (the dots 
represent the terms which do not contribute near the crit- 
ical point). The static critical exponents for the Gaussian 
model are 



Fourier transforming Eq. 



0,i 1,3,1,0 



(A7) 



5| gives 
a)rg(m^,q-) 



TgK^^ . (A8) 



Here, 'rn^_^, h^^^, Tq are the Fourier components of 
m{t, x), hit, x), r(|a;|), respectively. 

An interesting quantity in the critical phenomena is 
the response function Xu),q- 



Xu 



From Eq. 



8j), the response function is given by 



—iuj + {cq^ + a)Tq 



(A9) 



(AlO) 



The response function has a pole at a; = —i{cq^ + a)Tq, 
which implies that the relaxation time r„ behaves as 



icq' + a)Tq 



(All) 



In general, r^^p diverges near the critical point, which 
is the critical slowing down. This divergence is 
parametrized by a dynamic critical exponent z: 



Tq=o OC 



(A12) 



The dynamic universality classes are partly classified 
by the exponent z. We discuss two universality classes, 
model A and B. 

(i) Model A: In this case, Tq=o is a nonvanishing con- 
stant. Then, Eq. ((JTT|) gives Tq=o oc (T - Tj-i. The 
Gaussian model has ly ~ 1/2, so z 2. 

(ii) Model B: For model A, we have assumed that 
Fq^o 7^ Oi but this is not true for a system with a con- 
served charge, and the value of the exponent z changes 
as the consequence. If (?Ti^=o(i)) is a conserved charge. 



d{m^=o{t)) 
dt 







(A13) 



In order for this equation to be consistent with Eq. 
Fg^o ~ 0- The transport coefficient r(|a?|) has even par- 
ity, so Tq must be an even function of q. Thus, in the 
hydrodynamic limit, Tq (x q^ + 0{q^), and 

1 e 



To OC 



OC 



(A14) 



where the g-dcpendence is written as {^q); this scaling 
form can be justified from the dynamic scaling law. Then, 
z = 4. 

More generally, the critical slowing down can be shown 
from the dynamic scaling law. In general, the dynamic 
scaling relations for model A and B are 



2-7] (Model A) 
4-7] (Model B) 



(A15) 



Note that Eq. (|A15|) relates the dynamic critical expo- 
nent z to the static critical exponent 7]. For the Gaussian 
model, 77 = 0, so the relations (jAlSp reduce to the above 
obtained values for z. 

We discussed only model A and B, but the dynamic 
universality classes were classified by Hohenberg and 
Halperin [2J|, and they are known as model A, B, G, 
H, F, G, and J. These models are further classified by 
the values of dynamic critical exponents. 



B: Anomalous dimension for holographic 
superconductors 

As we saw in the text, the critical exponents of the 
order parameter are given hy v = 7/2, ?/ — 0, and z ~ 2, 
as long as the solution of Eq. ()4.4p can be expanded in a 
series of (tn, q). Then, most of the critical exponents are 
determined once the critical exponent 7 is obtained. In 
this Appendix, we show 7 = 1 by expanding the solution 
of Eq. (j4.4[) in a series of e^- Since 7 is a thermodynamic 
exponent, it is enough to consider a stationary homoge- 
neous perturbation ipo '■= V-'ro=q=o- 

Substituting ro = q = into Eq. p.ip . we obtain 



= 



pJLJLjL _ I 



du uP du 



2 2 
m 



H- 



2p 



f 



1+ K 



H 



Mu) , (Bl) 



where we used Eq. p. 31) . Also, a is defined in Eq. (|3.4p 
and is the control parameter for the holographic super- 
conductors.'' 

At the critical point (cr — CTc), there must exist a scalar 
hair solution of the equation 



d fc d ,22 

iP— /to 

du uP du 



~b7 



xP,{u) , (B2) 



^ The parameter k is related to cr by Eq. 1 13.51 1. 
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satisfying the boundary condition 



(B3) 



where the subscript "c" denotes the quantities at the crit- 
ical point. Note that the functional form of the electric 
potential 21 in the high-temperature phase must be valid 
even at the critical point by continuity. 

Let us solve Eq. (jBl|) in a cr-expansion. Defining e^- := 
1 — (j/(Tc, Sk := 1 — k/kc, and SipQ := ipo — ■0^, we obtain 



= 



d fc d ,2 2 

u'^— / m 

du uP du 



2p 

> -t J c 



1 



fc V 1 + -f^c 

O(e<,)-0,(u) + 0(Jk)V',(j/) + 0(e2 , e.^V^o) • (B4) 



Equation (jB4[) implies that jiJV'ol = O(ecr) since 5k oc 
(5cr cx e„ from Eq. p.Sp . Suppose that STpo{u) behaves as 



near the AdS boundary. Then, 



(B5) 



(B6a) 
(B6b) 



since V'o = '0c + ^i^o- Combining this with Eqs. (j4.14p 
and (|4.15p . we obtain the critical exponent of the ther- 
modynamic susceptibility 7 = 1: 



(B7) 
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